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Abstract. We have investigated the in- plane thermal conductivity Kab{T, H) of large 
single crystals of optimally oxygen-doped (Yi_2:,Sm:r)Ba2Cu307-5 (a;=0, 0.1, 0.2 and 
1.0) and YBa2(Cui_j^Znj,)307_5( y^O.0071) as functions of temperature and magnetic 
field (along the c axis). For comparison, the temperature dependence of Kab for as- 
grown crystals with the corresponding compositions are presented. The nonlinear field 
dependence of Kab for all crystals was observed at relatively low fields near a half of 
Tc- We make fits of the k{H) data to an electron contribution model, providing both 
the mean free path of quasiparticles Ho and the electronic thermal conductivity Kg, 
in the absence of field. The local lattice distortion due to the Sm substitution for 
Y suppresses both the phonon and electron contributions. On the other hand, the 
light Zn doping into the CuO 2 planes affects solely the electron component below Tc, 
resulting in a substantial decrease in £0 ■ 
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1. Introduction 

Since the discovery of high-Tc cuprate with its perovskite structure, a large number 
of measurements of thermal conductivity k have been reported up to date, giving 
crucial information about transport properties in the superconducting and normal 
states p. The thermal conductivity of solids is usually separated into the phonon 
and electronic components and He |2|- For high-T^ copper oxide superconductors 
YBa2Cu307_5(YBCO), Bi2Sr2CaCu208+5 and Tl2Sr2CaCu208+5 , a rapid enhancement 
in K,{T) is observed below Tc, with decreasing T, then k(T) reaches a maximum around 
a half of Tc and finally shows a monotonous decrease at low temperatures. Tewordt 
et al. ascribed the origin of the anomaly to the phonon contribution Kph on the basis 
of the classical theory of thermal transport in superconductors [21 • On the contrary, 
Yu et al. have proposed a new interpretation of the anomalous enhancement in k{T) 
as being responsible for the electron contribution Ke Following their proposal, a 
strong suppression in the quasiparticle scattering rate below Tc causes the n(T) anomaly, 
which is consistent with an observed peak in measurements of microwave conductivity 
in YBa2Cu307_5 |5^. In addition, the magnetic field effect on the thermal conductivity 
in high-Tc superconductors enables us to examine the scattering processes of thermal 
carriers due to magnetic vortex lines Hj. In ref. [Hj, the authors developed a two 
component method to analyze the field dependence of the thermal conductivity and 
concluded that only the electronic contribution is field sensitive. In a similar way, a 
procedure to check the linear field dependence of inverse electronic thermal conductivity 
was proposed by Pogorelov et al. ^9j. However, the origin of the anomalous peak in 
k{T) and its field dependence are still open questions. 

In this paper, we present measurements of the in-plane thermal conductivity Hab 
of Sm-substituted single crystals (Yi_j.,Sm2,)Ba2Cu307_5 (x=0, 0.1, 0.2 and 1.0) as 
functions of temperature and magnetic field. To our knowledge, no measurement of 
thermal conductivity has been performed on A-site substituted single crystals YBCO 
except for polycrystalline samples [TIH ITT] . Here, we give some comments on the 
substitution effect of Sm ion on the crystal growth of YBa2Cu307_5 . In a previous 
work, the substitution of Sm for Y led to a high growth rate for single crystal growth, 
which is about two times greater with a Sm content x=0.12 than that of YBCO |12j . 
Moreover, by controlling the Sm content, one obtained high superconductivity for Sm 
substituted YBCO. 

The K,{H) profile is fitted using the electron contribution model, to estimate the 
mean free path of quasi particles io near a half of Tc in the absence of field jHl IT8] . 
Compared with the linear field dependence of the inverse Kg ,the square-root fit to the 
K (if) data is applied . In addition, we discuss the impurity substitution effect on the 
thermal conduction of large single crystals YBa2Cu307_5. 
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2. Experiment 

A series of single crystalline (Yi_a;,Sma;)Ba2Cu307_5 (a;=0, 0.1, 0.2 and 1.0) and 
YBa2(Cui_j,Znj^)307_5 (?/=0.0071) were grown by the crystal pulling method. Sm 
contents range from 0, 0.1(0. 119), 0.2(0. 215) and 1.0. As-grown crystals with x=0, 0.1, 
0.2 and ?/(Zn) =0.0071 were annealed at 520°C in the oxygen atmosphere (for 1 to 2 
weeks)to achieve the optimal Tc with 7 — 6 =~ 6.9, as shown in Fig|21 SmBa2Cu307_5 
(SmBCO) crystal was oxidized at 300°C for 200h, giving a higher Tc property. The slight 
degradation in Tc for SmBCO probably arises from the partial substitution of Sm for 
Ba site rather than oxygen defects because we annealed it in the oxygen atmosphere for 
a very long time. The oxygen content for the as-grown crystal with x=0 was estimated 
to be 7 — 5 =~ 6.17 using the iodometric analysis on one piece cut from the same batch. 
Sample preparations are described by Yao et al.in [T21 CH US]- Sample dimensions 
are typically 4.5x3.1 mm^ in the a6-plane and 1.7 mm along the c-axis. The EPMA 
analysis showed that Sm ion is uniformly distributed in single crystals. The inductively 
coupled plasma (TCP) technique revealed that single crystalline samples used in our 
measurements are close to the nominal compositions. 

Measurements of the thermal conductivity, in the a6-plane and along the c-axis, 
were done with a conventional steady-state heat-flow method^H]- The magnetic field 
effect on the a6-plane thermal conductivity was examined in the application of fields 
up to 14 T along the c-axis at the Tsukuba Magnet Laboratory, the National Institute 
for Materials Science and at the High Field Laboratory for Superconducting Materials, 
Institute for Materials Research, Tohoku University. A temperature gradient across 
the sample was determined using a differential-type Chromel-constantan thermocouple 
after a calibration procedure. The magnetization measurement was carried out using 
the SQUID magnetometer. The temperature variation of the resistivity was measured 
on the identical samples using a four-probe technique. 

3. Results and discussion 

The in-plane thermal conductivity Kab for single crystals (Yi_2.,Sma;)Ba2Cu307_5 (x=0, 
0.1, 0.2 and 1.0) and YBa2(Cui_yZnj^)307_5 (y=0.0071) is shown as a function of 
temperature in FigO For comparison, the out-of-plane thermal conductivity Kcwith x=0 
and 1.0 is also depicted in the FigO For the identical crystals used for our experiment 
, measurements of the low-field magnetization and the in-plane resistivity are presented 
in Fig|21 to check the superconducting properties. The resistivity data for the lightly 
Zn doped sample will be shown in the inset of Fig.3(b). 

First of all, the Hab of (Yi_j;,Sm2,)Ba2Cu307_5 showed a rapid enhancement upon 
crossing with decreasing temperature, while for Kc no remarkable variation was 
detected near Tc. The substitution of Y by Sm strongly suppresses a maximum of 
Kab from 160 mW/cmK at x=0 down to 50 mW/cmK at x=0.2 through 85 mW/cmK 
at a:;=0.1. However, the pronounced peak of Kab below Tc is still observed even for the 
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x=0.2 crystal with its low thermal value. For the SmBCO crystal, Hab exhibits a low 
thermal conductivity accompanied by the broad maximum near 40K. This finding is 
closely related to a local lattice distortion due to the slight substitution of Sm for Ba as 
mentioned before. For free ions, the effective moment of Sm^"*" is 0.85 fiB and it is much 
smaller than the value of other rare earth ions (3.62 fiB for Nd^"*" and 7.96 hb for Gd'^"*"). 
If the phonon scattering due to magnetic ions with the smaller moment is assumed to be 
in the presence of the crystal field for SmBCO, it will contribute a suppression observed 
in Kab to some extent. 

Now, we comment on k in the normal state for the superconducting crystals. 
The thermal conductivity is separated into the electronic and phonon components. 
First, the electronic thermal conductivity in the normal state is estimated from 
the resistivity data using the Wiedemann- Franz (WF) law Kg = LqT / p with a Lorentz 
number Lq =2.45x 10~*Wf2/K^. For all superconducting crystals, the estimated value 
Kg is almost independent of T because of the T linear dependence of p, as shown in 
FigEfa). Accordingly, when ^ab — gives rise to the phonon component Kph{T), the T 
dependence of Kab{T) observed represents that of the phonon component in the normal 
state pTj . 

For comparison, let us show in Figl^l^b) the Kab data for the corresponding as- 
grown crystals of (Yi_2;,Sm2.)Ba2Cu307_5 (x=0 and 0.1) and YBa2(Cui_yZnj^)307_5 
(?/=0.0071) . In the inset of Fig. 3(b), the as-grown p data for the Zn-doped sample 
are presented and the electronic component is neglected from the WF law using 
its resistivity data. For other as-grown crystals, we also expect an insulating property 
or a weakly metallic one with high resistivity. We note that the as-grown crystals for 
our experiment were prepared at least under the similar crystal growth environment. 
Accordingly, we conclude that for all as-grown crystals used, k ~ Kp^. 

The Kab at x=0 follows a weakly T^^-like behavior, while for the Sm - substituted 
crystal it is almost independent on T. This finding is probably caused by the enhanced 
phonon scattering due to the Sm substitution. Here, At high temperatures, the T 
dependence of K.ab{T) of the as-grown crystal is expressed as Uph = l/(W^o + C(T), where 
Wq and a are thermal resistances due to defect scattering and phonon-phonon Umklapp 
scattering, respectively. Moreover, for the light Zn doping for the as-grown sample we 
observed a similar T dependence of k as that in the pure sample, indicating that a small 
amount of Zn-impurity does not affect the phonon thermal conduction. In contrast, 
the light Zn doping for the superconducting crystal strongly suppresses Tc down to 80K 
through the superconducting pair breaking, resulting in a reduced peak in Hab{T) in 
Figlla)^!!. 

We notice a common trend in the temperature dependence of phonon component 
between as-grown and superconducting samples. For example, i^abiT) of both as-grown 
and oxygen-annealed samples with x=0 shows a weakly T~^-like behavior at high-T. In 
addition, the oxygen content does not affect the T dependence of Kph strongly except for 
highly insulating crystal YBa2Cu306+5(5 < 0.11) [IHIIIII. Accordingly, let us assume 
that a minor modification of Hab{T) of as-grown crystals expresses the temperature 
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variation of the phonon component below for superconducting crystals if we ascribe 
the K peak to the electronic contribution. Thus, the original k data for as-grown crystals 
are scaled to coincide with Hab — i^e corresponding superconducting crystals at 

high temperatures (150K). The solid and dashed curves in FigQ represent the phonon 
contribution and , at the selected temperature of 40K, we obtain 43 % at a; = 0, 55 % at 
X = 0.1 and 71 % at ?/(Zn)=0.0071. (the corresponding electronic component is listed 
in Table 1) This rough estimation in is comparable with the value (~ 40% near Tc/2 
for YBa2Cu307_5 with 7 — 5 = 6.93 ) evaluated by Takenaka et al., [TJj. Separately, 
a primitive estimation on the phonon component of an untwinned crystal YBCO with 
Tc=90.5K gives a similar value jl]. 

Next, we show in the Figl^the magnetic field dependence of the in-plane thermal 
conductivity Kab{H) of single crystals (Yi_a;,Smj,)Ba2Cu307_5 (x=0, 0.1 and 1.0) and 
YBa2(Cui_j,Znj^)307_5 ( ?/=0.0071). The Kab{H) data are normalized by the zero field 
value. At low fields, the value of Kab exhibits a rapid drop around nearly Tc/2 and it 
then tends to saturate at high fields up to 14 T. The applied field causes a substantial 
decrease in k{H)/k{0) by ~ 30% at x=0, while k{H)/k{0) of Sm-substituted samples 
shows a smaller degradation from ~ 16% at x=0.1 to ~ 12% at x=1.0. In the mixed 
state of conventional superconductors, we beheve that thermal carriers (quasiparticles 
and phonons) are scattered by the magnetic vortices, resulting in a substantial decrease 
in K. 

Now, we try to fit the Kab{H) data to the electronic contribution model [Hlin]- The 
electronic thermal conductivity in magnetic fields, Ke{T,H) , is rewritten using both 
He(T, 0)and io as follows; 

..(T.i/) ^ ICMT,H) ^ ..(T.O)f^ ^ TTWmH) 

Here, the total mean free path of quasiparticles is expressed such as £(T,H)^^ = 
io{T)^^ +iy{H)~^, where io and are the mean free path of quasiparticles in zero field 
and in sufficiently strong fields, respectively. 

Accordingly, the field dependence of the thermal conductivity ,k{T,H), is fitted 
using the following formula 

k{T, H) = K^{T, H) + Kp,{T) = ^ + ^AT) (2) 

,where we set p(T) as cr£o/0o and ^o/^"\/0O) for ^^=1 a,nd 0.5, respectively {(po, the 
quantum flux). In the case of n = 1, a represents a scattering cross section per 
vortex cores in transport [20] • When the quasiparticle scattering from vortex cores is 
proportional to the number of the magnetic quantum flux, we get H)~^=aH/6n |21j. 
Thus, for the exponent n=l in eq. we take into account for the mean free path 
of quasiparticles iv{H) from the scattering due to vortex cores, in addition to io for 
scattering processes in the Meissner state. Using the a = 90 A in YBCO |2.1j, iv{H) 
is ~ 2.3 X lO^H^^ AT^^. On the other hand, in the case of n=0.5, assuming the 
quasiparticle scattering due to screening currents around vortex cores in a disordered 
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vortex array, the vortex contribution iy{H) is proportional to the average lattice constant 
of a vortex lattice a„ i.e., iv{H) = ka^ ~ k^j4>Q/ H, where is a constant independent 
of fields Pj. Using A;=0.5 in 1^{H) becomes ~ 2.2 x lO^if-Vs AT'^/^^ 

Our reliable fits of the n{H) data to Eq.© give the two parameters, Ke(T, 0)/k(=1- 
^^vh|^^ and p(cx ^o), at selected temperatures as summarized in Table 1. It is truly that 
our fitting procedures yield a better fit for the n = 1 than that for the n = 0.5 , but 
it is difficult to distinguish a if-linear or -\/ff-dependence in the denominator of k,{H) 
in Fig|3 In iY-linear fits, we notice the underestimation in electron component for all 
crystals if we accept the phonon component estimated from the as-grown data. For 
example, in the pure crystal of a;=0, the if-linear fit yields Kg ~ 32%, while in the \/H 
fit we get Kg ~ 50%. This tendency in Kg below Tc appears in previous works on the 
field dependences of k{H) for underdoped YBa2Cu306.63 and Tl2Ba2Cu06|Hl HI]- Such 
the underestimation in Ke will enable us to consider a significant modification of the 
free-electron Lorenz number in the WF law for the strongly electron correlated system 
such as high- Tc cuprates, giving the reduced electron thermal conductivity |23]. 

Finally, let us examine the fitting parameter p listed in Tabl^ From p cc io, 
we obtain the mean free path of quasiparticles io ~ 1.1 x 10^''m for n=l at the pure 
crystal. Here, we used the cross section of cr = 90 A in YBCO 21 . In a similar way, 
for the n=0.5, io ~ 7.8 x 10~^m with k=0.5 in [^. Both the former and latter fits 
give io ~1000A at 40K for pure x=0 crystal. From Tabl^ we guess that io at the 
Sm substituted crystal is about one-third as long as the value of io of the pure crystal, 
whether we use the if-linear fit or the If the normal state Ke(T) is proportional 

to •^o(^) at T > Tc , then the plots in Figinfa)give the ratio of io at x=0.1 to that 
at x=0 ~ 0.26 at higher T , This value is not far from the above results in the io ratio 
(~ 1/3) around a half of Tc. Moreover, for the lightly Zn-doped crystal, our fits to 
the data indicate the reduced mean free path of quasiparticles. A pair-breaking effect 
due to Zn-doping makes quasiparticles increase, giving the strongly electron-electron 
scattering enhanced below T^. This finding has a close relationship with the reduced 
peak observed in of YBa2(Cui_j,Znj^)307_5 with y=0.0071. 

4. Summary 

We have studied the in-plane thermal conductivity of large single crystals 
(Yi_^,Sm^)Ba2Cu307_5 (x=0, 0.1, 0.2 and 1.0) and YBa2(Cui_j,Znj,)307_5 ( y=0.0071) 
as functions of temperature and magnetic field. The nonlinear field profile of k,{H) 
observed around a maximum peak is discussed on the basis of the electron contribution 
model. We tried to make the fits of the fv{H) data to the if-linear or \^ dependence 
of the inverse electronic thermal conductivity in Eq.(2). The local distortion due to 
the Sm substitution for Y site suppresses both the phonon and electronic contributions, 
accompanied by the substantial decrease in io at 40K in comparison with that of the 
x=0 pure sample. The light Zn doping for Cu02 planes does not contribute the phonon 
conduction but affects solely the electron component below Tc, resulting in the reduced 
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mean free path of quasiparticles in the superconducting state. For SmBCO, a suppressed 
maximum in k,{T) arises from the partial substitution of Sm for Ba site ,which is closely 
related to the slight degradation in Tc. 

One of authors (X Yao) would like to thank NBRPCE(No.2006CB601003) 
and SSTC (No. 055207077) for financial support. We are grateful to F.Sato and 
Prof.S.Kambe for the iodometric analysis. 
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Figure 1. The in-plane tliermal conductivity Kab for single crystals 

(Yi_^,Sm^)Ba2Cu307-a (a;=0, 0.1, 0.2 and 1.0) and YBa2(Cui_yZny)307-a (y = 
0.0071) as a function of temperature. For comparison, the out-of-plane thermal 
conductivity Kcwith x=Q and 1.0 is also depicted. Solid and dashed curves in (a) 
represent the estimated phonon component Kph for x=0 and y = 0.0071, respectively 
(see to the text). In a similar way, a solid curve in (b) shows the phonon contribution 
Kph of the a;=0.1 crystal. 
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Figure 2. (a)normalized low-field magnetization M{T) under a zero-field 
cooling (ZFC) scan and (b) in-plane resistivity pab{T), for single crys- 
tals (Yi_a,,Sni3;)Ba2Cu307_5 {x=0, 0.1, 0.2 and 1.0). The M(T) data of 

YBa2(Cui_j,Znj,)307_5 {y = 0.0071) arc given in the inset of (a). The inset of (b) 
represents magnified plots for the Pab{T) data of the Sm-substituted samples. 
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Figure 3. (a) Normal-state electronic thermal conductivity estimated from 
the resistivity data using the Wiedemann-Pranz (WF) law k" = LqT/ p with a 
Lorentz number Lq =2.45x10~^WO/K^. (b) The Kab data for the corresponding as- 
grown crystals of (Yi_3;,Sm2;)Ba2Cu307_5 (.t=0 and 0.1) and YBa2(Cui_yZnj,)307_5 
(2/=0.0071). In the inset of (b), the Pab{T) data for the as-grown and oxidized samples 
with light Zn doping are presented. 
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Figure 4. The magnetic field dependence of tiie in-plane thermal conductivity 
Kab{H) of single crystals (Yi_3;,Sma;)Ba2Cu307_5 (a;=0, 0.1 and 1.0) and 
YBa2(Cui_yZnj,)307-(5 ( y=0.0071), at selected temperatures. The applied field H is 
parallel to the c-axis. Solid curves represent better fits of the n{H) data to the 7J-linear 
dependence of the inverse electronic thermal conductivity in Eq.(2). For comparison, 
the inset shows the curve fits to the \/^-dependence of the inverse Ke in Eq.(2). 
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Table 1. The fitting parameters Kc/k and p from tlie fit of Kab{H) data to eq.|21)for 
single crystals (Yi_^,Sm^)Ba2Cu307_5 (a;=0, 0.1 and 1.0) and YBa2(Cui_yZnj^)3 
Or-s (y^O.OOTl). For comparison, we show the electronic component near Tc/2, 
estimated from the k data of as-grown crystals described in the text. On the last 
column, the electronic components in the normal state are given using the WF law. 
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